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ABSTRACT 

Even since the first discovery of Poynting and Robertson, the radiation source has been treated 
as merely a point. Even in a very few studies where the size of the source has been taken 
into account, the treatment of the problem remained largely non-relativistic. In the present 
work, we address the issue of the finite size effects on the Poynting-Robertson effect in a fully 
relativistic manner for the first time. As a result, the emergence and the characteristic of the 
critical point/suspension orbit can be studied in a systematic and detailed manner. 
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1 Introduction 



Effects of the radiation field on a test particle orbiting a radiation source has been studied 
since the discovery of Poynting- Robertson effect [1] and [2] . We begin with the summary of the 
present status of research in the literature along this line. First, Guess [3] studied the effects of 
finite size and found that radiation drag (Poynting-Robertson effect) gets enhanced compared 
to the case of point-like radiation source. In his work he confined his calculation to first order 
in v/c. Second, Carrol [4] found that the finite size effect of the radiation source renders 
the radiation drag force greater than its counterpart for a point-like source. His treatment, 
however, was limited to the context of special relativity. Third, Abramowicz, Ellis, and Lanza 
[5] constructed the general relativistic radiation stress-energy tensor describing the radiation 
field from the non-rotating radiation source with finite radius size for the first time and solved 
the equation of motion to discover the existence and the location of a critical point at which all 
segments of accretion flow are captured. As their concern was to describe jets and outflows, they 
limited the particle's motion to one-dimensional radial motion alone. Fourth, Miller and Lamb 
[6] and [7] considered and solved the azimuthal as well as radial component of the equation 
of motion. They, however, failed to notice the emergence of and study the nature of critical 
point as they considered the set-up in which the luminosity of the radiation source is well below 
the Eddington critical value. Interestingly enough, lastly, Bini, Jantzen, and Stella [8] studied 
general relativistic version of the Poynting-Robertson effect by considering and solving all 3- 
components (time, radial, and azimuthal) of the equation of motion. Their set-up, however, 
is not identical to that in our present work in that they employed the radiation stress-energy 
tensor describing the radiation field from a point source. They ignored the finite size of the 
radiation source. 

In the present work, in order to address general relativistic version of finite size effects on 
the Poynting-Robertson effect in a rigorous and complete manner, we employed the radiation 
stress-energy tensor first constructed by Abramowicz, Ellis, and Lanza [5] to explore the detailed 
impact of the radiation forces on the test particle's trajectories. 

In section 2 we describe the derivation of the equation of motion of the test particle in the 
Schwarzschild space-time background outside a non-rotating, isotropic radiating spherical star 
with finite size. In section 3, we present the results of the numerical integrations along with 
the our analytical results. Finally in section 4, we summarize and discuss our results. 

2 Equations of motion 

To derive the equations of motion that govern the orbital motion of the test particle in the radia- 
tion field from the spherical star emitting the radiation isotropically, we start with Schwarzschild 
spacetime, 

ds 2 = g^dx^dx" 
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= -(1 - 2M/r)dt 2 + (1 - 2M/r)- 1 dr 2 + r 2 [dQ 2 + sin 2 ^0 2 ). (1) 

We work in the geometric units, where G = 1 = c (G is the gravitational constant and c is the 
speed of light). The equations of motion is given by 

a a = — (2) 
m 

where f a denotes a non-gravitational force (for example, radiation force) exerted by the radi- 
ation (or luminosity) on the test particle, m is the rest mass of the test particle, and 

dU a 

a a = -^r + T^lFir (3) 

is the acceleration, where U a is the four-velocity of the particle, T^ v = \g al3 {gi3^,v + 9pu,^ — 9nv,p) 
is the Affine connection, and comma (,) denotes partial derivatives. 

If we assume that the cross-section a of the test particle is independent of energy (frequency) 
and direction of the photon, the non-gravitational force f a due to scattering of photon from 
the star is given by [9] , 

f a = oF a , (4) 

where F a shown below is the quantity transformed from the radiation energy flux measured 
in the particle's rest frame (which is comoving with the particle) using the orthonormal tetrad 
e? associated with the particle's rest frame as follows, 

r ~ i 1 CO 

= -h«T^U a , (5) 

where hp = 8p + U a Up is the projection tensor, and is the radiation stress-energy tensor 
in the lab frame. 

According to [5], the components of the radiation stress tensor T a/3 measured in LNRF 
(Locally No n- Rotating Frame; [10, 11]) can be written as, 



T tt = 27r/(r)(l -cos a) 


(6) 


T if = 7il(r)sm 2 a 
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T H = -7r/(r)(l -cos 3 a) 


(8) 


T ee = ^/(r)(cos 3 a -3 cos a + 2) 


(9) 


T** = -J(r)(cos 3 « -3cosa + 2), 
o 


(10) 



where a is an apparent viewing angle of the star seen by an observer (FIDO; fiducial observer) 
at rest in the LNRF and is given by sin a = (jj^j ( II^m/r ) ( see ^ or ^ ne ra dius of the star 



R > 3M, and I(r) is the frequency-integrated specific intensity at the radial position r and can 
be expressed as (see Appendix A in the [7]) 

(1 - 2M/R) mM ( L°°\ 
{ > (l-2M/r)%,ff UfJ' 1 > 

where L°° is the luminosity at infinity and L^ dd is the Eddington critical luminosity which is 
given by L^ dd = AnmM/a. 

By transforming T^ v using the tetrad associated with the LNRF, we can obtain the 
components of the radiation stress tensor T a/3 , 

T af) = efeZTi*, (12) 

where the tetrad associated with LNRF in the Schwarzschild coordinates is given by, 

e 5 = (1 - 2M/r) 1 ' 2 dt 1 

e 1 = (l-2M/ry 1/2 dr, (13) 

e 2 = rd6, 

e 3 = rsin9d(f). 

Finally, the radiation stress-energy tensor describing the radiation field from the radiation 
source of finite size (R > 3M) can be rewritten as, 

*--(7)^®( i ~HS 

-e)(")(-TfO 

T „ = f m \ 2(l + cosa + cos 2 a) /M\ ( L°° \ ^ 



a J 3(1 + cos a) \r 2 J\L°° 



EM, 



T »» i m ) (2 -cos a -cos 2 a) /M\ / 2M- ^ 



a/ 3(1 + cos a) Vr 4 / V r ) \L^ dd / 
H fm\ (2 -cos a -cos 2 a) f M w 2MX' 1 / \ 

W 3(l + cosa) vr 4 sin 2 #A r) \L% d J ' 1 J 

When a radiation source is point-like, the apparent viewing angle (a) measured by FIDO 
is equal to zero, thus cos a = 1. Inserting this cos a = 1 into above equations (14) through 
(18) gives the radiation stress-energy tensor for point-like radiation source, which is exactly 
identical to the radiation stress-energy tensor employed by Bini, Jantzen, and Stella [8]. 

When the radius (R) of the radiation source is smaller than the photon sphere (3M), i,e., 
R < 3M, the apparent viewing angle is given by, 

"M\ / 2M\ 1/2 



sin a 



r /M\ / 2M\ 1/Z . . 

= 3V3(-)(l--) . (19) 
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Therefore, the radiation stress-energy tensor describing the radiation field from the radiation 
source which has a radius smaller than the photon sphere have no dependence of the source's 
radius R. Thus, such radiation sources can be treated as a radiation source with R = 3M(the 
radius of the photon sphere). 

We can restrict the orbital motion of the particle to the equatorial plane (9 — f , Uq — 0), 
without loss of generality owing to the spherical symmetry of non-rotating spherical star, thus 
the equation of motion finally can be decomposed into each components as follows; 



dUt 
dr 



(£)[l-2(l-^)~'u? 

(8 + 2 cos a + 2 cos 2 a) /M\ 

3(1 + cos a) W 
(8 — cos a — cos 2 a) 



r" ) 



3(1 + cos a) 



T OO 

> ^Edd , 



T OO 

^Edd, 

2M^ 



U t U 2 r 



T OO 

^Edd, 



(20) 
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cos 2 a) 
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(21) 



U r 



dUf> 
dr 



— cos a — cos a 
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(22) 



We are ready to envisage the features of the solution to the equation of motion. First, 
equation (21), radial component of equation of motion, includes four lines: the three terms in 
the first line is not relevant to the radiation, as overall sign of the second line is positive definite, 
this line can be identified with the outward radiation pressure gradient and it should be noted 
that the second line has no dependence of the finite size of the radiation source. Both third line 
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and fourth line can be regarded as being responsible for the radial radiation drag force as their 
overall sign are negative definite. As the third line depends on the radial velocity component 
alone, if we take into account only radial motion, the radial radiation drag is attributable to the 
third line. As we can know from the finite size dependent part ( ^-^^^~^ ~ ) of the third 
line, the radial radiation drag at the surface of the radiation source (cos a = 0) gets 4/3 times 
larger than that of the point-like source (cos a = 1), which is identical with the results in the 
study of Guess [3] and Carrol [4]. However, their study failed to mention about the our fourth 
line radiation drag which depends on both the radial velocity and the azimuthal velocity. The 
radial radiation drag in the fourth line gets 8/3 time large compared to the point-like source. 
Second, equation (22), azimuthal component of equation of motion, is composed of three lines: 
the first line depends on both the radial and the azimuthal velocity and its overall sign relies 
on the direction of the radial motion of the particle. When the particle is in inward motion, 
as its overall sign is negative definite, the first line can be regarded as the azimuthal radiation 
drag, but if the motion of the particle is changed to outward one, as its overall sign is changed 
to positive definite, it can be identified with the radiation counter-drag. As the second line has 
a dependence of the azimuthal velocity only and its overall sign is negative definite, this second 
line can be regarded as the azimuthal radiation drag exerted on the particle in circular motion. 
The finite size dependent part in the second line indicates that the azimuthal radiation drag 
at the surface of the radiation source gets 8/3 times larger than that of the point-like source, 
which also is identical with the results in the study of Guess [3] and Carrol [4]. Finally, the 
third line depends on both the radial and the azimuthal velocity and its overall sign is negative 
definite, thus this third line can be also identified with the radiation drag. Interestingly, the 
finite size of the radiation source leads to smaller radiation drag than a point-like source. 



3 Results of numerical integration 

For the sake of convenience, we introduce luminosity parameter L = -p^-. 

^Edd 

Figure 1 shows the ways the test particle arrives and enters onto the dotted circle. The 
test particle arriving at the dotted circle seems to be suspended and not to move any more 
there, thus hereinafter we will refer to the dotted circle as "suspension orbit". The starting 
points of the particle are the same as (6M, 0) in Cartesian coordinates (X,Y). Initial azimuthal 
velocities are fixed to be 0.2c and its direction is counterclockwise. Also, the luminosity of the 
radiation sources are the same as L = 0.8. Radial distance of the "suspension orbit" is roughly 
r so ~ 5.56M. Solid line denotes the trajectory of the test particle when the radiation source 
is point-like, and dotted, dashed, and long dashed lines denote the trajectories of the particle 
when the radiation sources have the radii of 3M, 4M, and 5M, respectively. As it is mentioned 
above section II, it should be noted that the radiation source with radius smaller than 3M 
(photon sphere) can be treated as one with radius of 3M. 

As shown in Figure 1, as the size (radius) of the radiation source gets larger, the "suspension 
orbit" stays the same (location) but the test particle arrives at the "suspension orbit" more and 
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X-axis (M) 



Figure 1: shows the ways the test particle arrives and enters onto the suspension orbit. 4- 
trajectories provided exhibit as the radiation source gets larger and hence it's surface gets 
wider, the Poynting-Robertson effect gets more and more manifest reducing the test particle's 
azimuthal speed more and more remarkably. The luminosity parameter of the radiation sources 
are the same as L = 0.8 and dotted circle denotes "suspension orbit" . 




-0.015 r 

5.6 5.8 6.0 6.2 6.4 
Radial distance (M) 



Figure 2: demonstrates the radial velocity profiles of the 4-trajectories of the test particle given 
in Figure 1. 
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Figure 3: demonstrates the azimuthal velocity profiles of the 4-trajectories of the test particle 
given in Figure 1. 

more quickly. Unlike the case of the point source where the test particle hits the "suspension 
orbit" at some finite angle, when the source has finite size, the test particle always strikes the 
"suspension orbit" nearly perpendicularly. 

Figures 2 and 3 demonstrate the radial and azimuthal velocity profiles of the 4-trajectories 
of the test particle given in Figure 1, respectively. As shown in Figures 2 and 3, around the 
"suspension orbit", as the size of the radiation source gets larger, the radial and azimuthal 
velocities are more and more quickly reduced and finally vanish at the "suspension orbit". 
Thus, from Figures 1 through 3 we can find out that the radial and azimuthal velocities of the 
test particle at the "suspension orbit" are, respectively, equal to zero and the rate of change of 
the radial and azimuthal velocities also vanish, respectively, hence the test particle comes to a 
complete rest and does not move any more there. In other word, the test particle there satisfies 
the following conditions, 



U r = f^ = 

dU t dU r _ dUj 
dr dr dr 



For U r — — 0, the normalization condition g^U^U^ = — 1 yields U t — — (1 — 2M/r) 1 ' 2 . 
Equations (20) and (22) satisfy the above conditions, and Plugging this U t — — (1 — 2M/V) 1 / 2 



together with = into equation (21) gives, 
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(M\ ( 2M\- 1/2 (M\( 2M\~ 1 (L co \ . . 

» = -b)(i-— ) + b)( 1 --) (z=J- < 23 > 



The solution of Equation (23) gives the radial distance r so of the "suspension orbit" , 



2M 



1 - 



J oo 
, ^Edd , 



(24) 



which is equivalent to the radial position of the critical point in Abramowicz et al.[5] and Bini 
et al.[8], and is determined by ( ) alone. 

V ^Edd J 

Since first term on the right-handed side in equation (23) denotes gravitational inward force 
exerted on the test particle at complete rest in the "suspension orbit" and second term refers to 
the outward radiation pressure there, equation (23) means that the gravitational inward force 
at the "suspension orbit" exactly balance the outward radiation pressure there. 

It is interesting that among the entries/components of the radiation stress-energy tensor, 
the (t,t)-component, which is the radiation energy density and the (r,r)-component, which is the 
radiation pressure, exerts only on the moving particle, while the only (t,r)-component, which is 
the pressure gradient (force) exerts on both the stationary particle and the moving particle. As 
we already mentioned above, the (t,t) and (r,r) components of our choice of the stress-energy 
tensor are different from those of Bini et al.[8] whereas the (t,r) component of ours is exactly 
same as that of Bini et al.[8]. 

As a result, although the trajectory of the test particle before it reaches the "suspension 
orbit" in our study is manifestly different from that in the study of Bini et al.[8], once the test 
particle arrives at the "suspension orbit" , the outward radiation pressure gradient that balance 
the inward gravitational force in our study is the same as that in the study of Bini et al.[8] and 
this is precisely why our radial distance of the "suspension orbit" and theirs turn out to be the 
same. 



4 Concluding remarks 

Although since the remarkable discovery first made by Poynting and Robertson, a number 
of study on this Poynting-Robertson effect have been performed, in almost all the works the 
luminous radiation source has been treated as being point-like. This could be due to their 
particular aims of investigations in which they just tried to eliminate unnecessary complexity 
which has little to do with essence of their studies. Largely, however, such simplification or 
reduction of the set-up can be attributed to the technical and computational difficulties. In 
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order to render the set-up for the careful and detailed investigation of the Poynting-Robertson 
effect to be more realistic and practical, obviously one has to take into account the finite 
size of the radiating source. Even in a very few works in the literature addressing this issue 
of finite size, the treatment of the issue has been very limited. For instance, Abramowicz 
et al. [5] attempted general relativistic and analytical approach, but they considered radial 
equation of motion alone. Next, Miller and Lamb [6,7] performed general relativistic and 
numerical approach, but they considered angular equation of motion alone. To summarize, 
even these few works addressing the issue of finite size effect for the Poynting-Robertson effect 
have remained essentially incomplete. Along this line, our present work can be regarded as 
comprehensive treatment of the finite size effects on the Poynting-Robertson effect in a fully 
general relativistic context for the first time ever. As a result of this comprehensive treatment of 
the problem, we were successful in realizing the existence of critical point/suspension orbit and 
in understanding the detailed characteristics of this suspension orbit. It is our hope, therefore, 
that the present study of ours could solve as a ground work on which further serious studies 
addressing the issue of finite size effects on the Poynting-Robertson effect can be successfully 
carried out. 
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